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The mechanism of the magneto-Coulomb oscillation in ferromagnetic single electron transis- 
tors (SET's) is theoretically considered. Variations in the chemical potentials of the conduction 
electrons in the ferromagnetic island electrode and the ferromagnetic lead electrodes in magnetic 
helds cause changes in the free energy of the island electrode of the SET. This is a plausible 
origin of the conductance oscillation of the SET in sweeping an applied magnetic field. 
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§1. Introduction 

There have been intensive studies on the single electron charging effects^ & last several years. In 
experimental works on metallic devices, however, the materials investigated have been restricted to 
a few kinds of metals such as Al, for example. Combination of materials of different characteristic 
properties, such as normal-conducting metals, superconductors, ferromagnetic metals, could add 
more fertile aspects to the physics of the single electron devices. 

Recently the present authors have found several novel phenomena in ferromagnetic single electron 
devices.&B* One of them is the magneto-Coulomb oscillation in a ferromagnetic single electron 
transistor (SET). lit The ferromagnetic SET in ref. 4 is composed of a Co island electrode and Ni 
lead electrodes. A conductance oscillation similar to the Coulomb oscillation^ was observed in 
sweeping an applied magnetic field. In this paper we discuss the mechanism of this phenomenon 
with supplementing experimental results on a Ni/Co/Ni ferromagnetic SET. 

The essential origin of the magneto-Coulomb oscillation in the ferromagnetic SET is thought of 
to be noticeable change in the chemical potential of an isolated ferromagnet in magnetic fields. This 
causes different influences on the operation of the SET depending on whether its island electrode is 
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composed of a ferromagnetic metal or the lead electrode is ferromagnetic. We describe the magneto- 
Coulomb oscillation based on the linear response theory of the Coulomb oscillation!' @* taking into 
consideration the behavior of the ferromagnetic component of the SET in magnetic fields. Bearing 
the experimental situation in ref. 4 in mind, we consider the case in the classical regime where 
average spacing of the one-electron energy levels (Ae) in the island electrode is much smaller than 
the thermal energy of the electrons in the system (Ae <C k^T) and the energy distribution of 
the electrons in the island electrode is well described by the Fermi-Dirac distribution.!) We also 
assume that the tunnel resistance of the tunnel junctions are sufficiently larger than the quantum 
resistance, Rq(= h/e 2 ), and the orthodox theory of the single-electron tunneling^ is applied to the 
case without considering the quantum charge fluctuations or the renormalization of the charging 
energy. 

The paper is organized as follows: In §2 magnetic-field-induced variations in the chemical po- 
tential and in the work function of the ferromagnetic island electrode are examined. Here we 
regard the island electrode as an isolated small grain of a ferromagnetic metal. In §3, effect of 
the magnetic field for the ferromagnetic lead electrode and its influence on the operation of the 
SET are considered. Based on the consideration in these sections, the thermodynamic potential 
of the ferromagnetic island electrode is constructed in §4. This quantity dominates the charac- 
ter of the transport through the SET in the linear response regime.B& In §5, we examine the 
magneto-Coulomb oscillation in the ferromagnetic SET in some detail adopting a simple model for 
the ferromagnet. More realistic consideration on the ferromagnet and corrections to the results in 
§2-5 are presented in §6. Experimental results on the magneto-Coulomb oscillation in the Ni/Co/Ni 
SET which supplement the result in ref. 4 are presented in §7 and are discussed in §8. Possible 
applications of this phenomenon is discussed in §9, and the conclusion is described in §10. 

§2. Change in the Work Function of the Ferromagnetic Island Electrode 

The behavior of the ferromagnetic metal in a magnetic field has many aspects, and it has not been 
sufficiently understood until now. Therefore, in order to describe the essence of the mechanism of 
the magneto-Coulomb oscillation in the ferromagnetic SET, we adopt for a moment the simplest 
model for a ferromagnetic metal, that is, a degenerate free electron system that has spin-dependent 
density of states, and ignore the exchange-correlation effect for excitations. We also ignore effect 
of the orbital angular moments of the electrons. These contributions will be discussed in § 6. 
Moreover, for the sake of simplicity, we assume, in the following discussions, that the ferromagnet 
under consideration has a single domain. The effect of the multi-domain structure of the actual 
ferromagnet will be discussed later. 

Suppose a ferromagnet in an applied magnetic field, H, with the direction of its magnetization 
parallel to H. Because of the Zeeman effect, the energy band of the majority-spin electrons shifts to 
the lower energy, in comparison with the case in zero magnetic field, by 5(H) = gfisH/2, where g 
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is the electronic g-factor and /xb is the Bohr magneton. On the other hand, the energy band of the 
minority-spin electrons shifts to the higher energy by S(H). In equilibrium, re-population between 
the majority-spin electrons and the minority-spin electrons occurs through spin- flip scatterings 
(Fig. 1). Now we consider a ferromagnetic island electrode with electrons of a fixed number N. N 



Fig. 1. Schematic energy band diagram of the majority(+)- and the minority(— )-spin conduction electrons in a 
ferromagnetic metal in zero magnetic field (broken line) and in a finite magnetic field H (solid line). 2 A denotes 
the exchange splitting of the two subbands in zero magnetic field. 8(H) indicates the magnitude of the shift of 
the bands in H. ("(0) and C(^) respectively indicate the chemical potentials of the conduction electrons in zero 
magnetic field and in H. W(0) and W(H) means the work functions in zero magnetic field and in H, respectively. 
—£ c is the midpoint between the bottoms of the two conduction subbands measured from the vacuum level. 

is a sum of the number of the majority-spin electrons, N + , and that of the minority-spin electrons, 
The equivalence of the total number of electrons both in zero magnetic field and in magnetic 
field H is expressed as follows: 



Here we take, as the origin of the energy of an electron, the midpoint between the bottoms of the 
conduction subbands of the majority- and the minority-spin electrons, and name the difference of 
the energy between this point and the vacuum level of an electron S c (Fig. 1). Throughout this 
paper, the vacuum level is defined at infinity where magnetic field and the electrostatic potential 
is considered to be zero.H' In the above expression p±(e) indicate the densities of states of the 
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majority(+)- and the minority(— )-spin electrons in zero magnetic field, and f(x) is the Fermi- 
Dirac distribution function: 

= it^t- (2 - 2) 

The symbol A denotes a half of the magnitude of the exchange splitting at the bottoms of the 
conduction subbands and C{H) is the chemical potential of the conduction electrons in a magnetic 
field H%} The work function in a magnetic field H is expressed as W(H) = £ c — C(H)- 

Equation (2.1) is rewritten to a balance equation for the re-population between the majority- 
and the minority-spin electrons in magnetic fields: 

oo 

p+(e)[f(e - C(0)) - f(e - 5(H) - C(H))]de 

-A 

/oo 
p-(e)[/(e - C(0)) - /(€ + 5(H) - ((H))]de = 0. (2.3) 
,-A 

In ferromagnets such as 3d transition metals, a condition, k-gT <C C~ ^> is satisfied at temperatures 
below room temperature, and the lower limits of the above integrals can be replaced by 0. In not 
extremely strong magnetic fields, a condition, 5(H) /£ <C 1, is satisfied and we expand the Fermi- 
Dirac function as the power of 5(H) /((0) and A((H)/((0) (A((H) = ((H) - ((0)). Up to the first 
order in these quantities, eq. (2.3) becomes 

ro ° P+(e){S(H) + A((H)} - p-(e){5{H) - A((H)} 

l + cosh[(e-C(0))A B T] ' l ' J 

This equation determines the change in the chemical potential A£(H), and, therefore, that in the 
work function of the ferromagnet in magnetic fields. If the densities of states of the conduction 
electrons over an energy region of the order of k-gT around Q is constant, eq. (2.4) gives a simple 
relation: 

( P+ - p_)6(H) + ( P+ + P -)AC(H) = 0, (2.5) 
where p± = /?±(C(0))- This gives the shift of the chemical potential: 

A((H) = -Vgp B H/2, (2.6) 

where 

V = P -±^P-. (2.7) 

P+ + P- 

This quantity V denotes the spin polarization of the electronic densities of states at the Fermi 
energy of the ferromagnetic metal under consideration. It should be noted that this quantity is not 
the same one as measured in the spin-polarized tunneling experiment. The latter is weighted by 
the tunneling probability of each conduction band at the Fermi energy.0'0lll) 

§3. Effect of the Magnetic Field for the Ferromagnetic Leads 

The chemical potential of the ferromagnetic lead electrodes also changes in magnetic fields, which 
causes another effect on the operation of the SET. In real experimental situations, the ferromagnetic 
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leads are connected to non-magnetic metallic leads, which are connected to electronics, and the 
whole system is finally anchored to the ground somewhere in the circuit. In order to make the 
discussion clear, we assume a system depicted in Fig. 2 in the following discussions. 



Fig. 2. The ferromagnetic SET. A and B denote a ferromagnetic island electrode and ferromagnetic lead electrodes, 
respectively. The ferromagnetic metals of A and B is not necessarily the same. C denotes a non-magnetic metallic 
lead. C and C g respectively indicate the capacitance of each tunnel junction and that of the gate capacitor. V g 
is the electrostatic-potential-difference source for the gate and it also means the electrostatic-potential difference 
itself in the text. 

In general, there exists some difference in the electrostatic potential between two kinds of metals 
in contact, because the work functions of the two metals before contact are different in most 
of cases. Once they are in contact, electrons are transferred from the metal of the lower work 
function to the other, and an electric dipole layer is generated at their interface up to the depths 
around their screening lengths. The electric dipole layer makes their electrostatic potentials shift, 
and, therefore, changes the electrostatic-potential energies of electrons in the two metals so as 
to equate their electrochemical potentials (and the work functions) of both metals. Therefore 
electrostatic potentials of the ground, the non-magnetic leads and the ferromagnetic leads in Fig. 
2 are considered to be different in general, although the electrochemical potential, u, and the work 
function, W, are kept constant over the combined system in equilibrium. 

In magnetic fields, as describe in the previous section, the chemical potential of the ferromagnetic 
lead electrode changes by A((H) = —VgfJ,BH/2. Assume V is negative, for example, in the 
ferromagnet under consideration. Then AQ(H) becomes positive. In this case electrons flow from 
the ferromagnetic leads to the non- magnetic leads through their interfaces, the charges of the 
electric dipole layers at the interfaces change, the electrostatic potential of the ferromagnetic leads 
is raised, and the electrostatic-potential energy of conduction electrons is lowered to compensate 



5 



the change in its chemical potential. As a result, the contact potential difference, V c , changes by 

AV C (H) = A((H)/e 

= -Vgfi B H/2e, (3.1) 

and the electrostatic potential of the ferromagnetic leads is raised by AV C against the non-magnetic 
leads. The work function W of the combined system, however, remains constant irrespective of 
the magnitude of the applied magnetic field. It is because the nonmagnetic lead is connected to 
the ground, the difference of the electrostatic potential between the nonmagnetic lead and the 
ground is fixed, and, therefore, the electrochemical potential u measured from the vacuum level at 
infinity is fixed over the whole combined lead independent of the applied magnetic field. The overall 
change in the electron system of the combined lead caused by the magnetic field is a decrease in the 
electrostatic-potential energy of electrons in the ferromagnetic leads by — eAV c {H) (and, of course, 
the increase in the chemical potential by A((H)). 

The influence of this effect on the SET device is the change in the electrostatic potential of the 
island electrode against the non-magnetic lead electrodes through its capacitive coupling with the 
adjoining ferromagnetic leads. For example, in the situation of Fig. 2, the change in the electrostatic 
potential of the island electrode becomes AV C (H) ■ 2C/{2C +C g ), where C is the capacitance of each 
tunnel junction and C g is the capacitance of the gate capacitor. 

We should note the mechanism is effective even if all the leads except the gate electrode are 
ferromagnetic. This is because the spatial variation of the magnetic field causes the gradient of 
electrostatic potential in a ferromagnetic lead. The electrostatic potential of the parts of the lead 
electrodes, which are adjacent to the island electrode via the tunnel barriers, against the ground 
or against the gate electrode is determined by the magnetic field at the junctions irrespective of 
the kind of metals in between. If the gate electrode is ferromagnetic as well, the variation of its 
electrostatic potential due to the magnetic field should be included in the consideration. 

§4. The Thermodynamic Potential for the Island Electrode 

It has been indicated that conductance of SET-type devices at small bias voltages (eVbias *C k^T) 
are predominated by a distribution function, P cq (N), of the electron number in the island electrode 
in equilibrium MH* The distribution function is expressed by using a Gibbs thermodynamic potential 
Q(N) of the island electrode as§'& 

Now we consider the Gibbs thermodynamic potential of the island electrode of the ferromagnetic 
SET in a magnetic field. For the sake of clarity, we consider especially the case of ferromagnetic 
SET in Fig. 2; the system is composed of an island electrode of ferromagnetic metal A, adjoining 
lead electrodes of ferromagnetic metal B, other electric circuit made of non-magnetic metal C and 
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a source of electrostatic-potential difference for the gate. We assume the circuit is connected to the 
ground at the non-magnetic lead C and is placed in a space with no background gradient of the 
electrostatic potential. 

First we examine the electrochemical potentials of the majority- and the minority-spin conduction 
electrons in the island electrode A. The chemical potentials of the majority- and the minority-spin 
electrons, £±(H), measured from the bottom of each conduction band, are expressed as 

C ± (H) = C(H)±A±5(H). (4.2) 

If the shift of the electron bands relative to the vacuum level due to the charging effect is not 
considered, the energies of the bottoms of the conduction subbands, £±, measured from the vacuum 
level, are expressed as 

S± = -£°tAtS(H), (4.3) 

where £® means the energy distance between the vacuum level and the midpoint between the 
bottoms of the conduction bands without the charging effect. From eqs. (4.2) and (4.3), the 
electrochemical potentials of the majority- and the minority-spin electrons, uj±(H), measured from 
the vacuum level, become 

u±{H) =(±(H)+£± 
= ((H) - £° c 

= u A {H). (4.4) 

The electrochemical potential uja(H) is equivalent to — Wa(H), where Wa(H) is the work function 
of the ferromagnet A in magnetic field H. Therefore Gibbs free energy of the conduction electrons 
without consideration of the charging effect becomes 

N + u + {H) + AT_w_(F) = Nlu a (H). (4.5) 

It should be noted here that now energy of an electron is measured from the vacuum level at infinity. 

Next we take the charging effect into account. If the self capacitance €■£ of the island electrode is 
small, a change in the electron number from that of the electrostatic neutrality causes a considerable 
change in the electrostatic potential of the island electrode relative to the lead electrodes and 
electrostatic energy is accumulated in the system. Electrostatic-potential energy of conduction 
electrons is also changeable homogeneously by external sources. We label the electrostatic energy 
of the island as U, which includes both components due to the charging effect and external sources. 
In the case of Fig. 2 the external source is the electrostatic-potential-difference source for the 
gate,© and U is a function of N, V g and H. Then the Gibbs free energy of the island electrode 
becomes Noj a (H) + U(N,V g , H), and the Gibbs thermodynamic potential, Q(N, V g , H), in this 
case is expressed as 

n(N, V g , H) = Nloa(H) + U(N, V g , H) - Noj bc , (4.6) 
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where wbc denotes the electrochemical potential of the combined leads B and C measured from 
the vacuum level. Magnitude of wbc is equivalent with the work function of the combined system 
B and C, Wbc, which is independent of the applied field, as discussed in §3. 

In the present discussion, we assumed wa (k>±) has no N dependence, which is reasonable because, 
although the metallic island electrode is of the mesoscopic size, it still contains a large number of 
electrons. It results in a dense density of states at the Fermi energy. For example, in Ni of 
dimensions of 100nmx2/imx20nm, the density of states at the Fermi energy is about 9xl0 9 states 
Ryd -1 (see Table I). Thus a change in N of unity causes a change in the chemical potential of the 
order of 1(T 9 eV (~ l(r 5 K), which is negligibly small compared with the magnitude of a typical 
charging energy (~ IK). 

For the sake of clarity and completeness of the discussions, we hereafter consider the polarity 
of the magnetic field explicitly, and introduce spatial polarity into the spin polarization of the 
electronic densities of states. We take the direction of positive magnetic field as z direction and 
define a spin polarization with spatial polarity, V' , as 

+V ( magnetization is parallel to z ) 
V = I (4.7) 
I —V (magnetization is antiparallel to z). 

According to this notation, V'H becomes V\H\ if the magnetization and the applied magnetic 
field is parallel and becomes — V\H\ if they are antiparallel irrespective of whether H is positive or 
negative. We call V' simply as the 'spin polarization' hereafter. 
Using this notation, u A {H) is expressed, according to the consideration in §2, as follows: 

u A (H) = lo a (0) - V' A g A mH/2, (4.8) 

where V A denotes the spin polarization of the ferromagnet A in zero magnetic field. 

0) 

The elec- 
trostatic energy U(N,V S ,H) is expressed as follows: 

(N - N n ) 2 e 2 

U(N, V s , H) = ± ^ Ne0 ext (V g , H), (4.9) 

where Nq denotes the number of electrons when the island electrode is electrostatically neutral, 
Cs = 2C + C g and 

C 2C 

4>e*{V g ,H) = ^Vg + -K(ff). (4.10) 

Here V g is the electrostatic potential of the gate electrode relative to the non-magnetic lead C due 
to the electrostatic-potential-difference source for the gate and V C (H) denotes the contact potential 
of the ferromagnetic leads B relative to the electrostatic potential of the non-magnetic lead C in 
magnetic field H. According to the discussion in §3, V C (H) is expressed as follows: 

V C (H) = V c (0) - V^g Bf i B H/2e, (4.11) 
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where V c (0) is the contact potential of the ferromagnet B against the non-magnetic metal C in zero 
magnetic field, and V B denotes the spin polarization of the ferromagnet B in zero magnetic field. 
After rearrangement of all the entities in 0(iV, V g ,H), we obtain the next expression: 

,2 



Q(N,V g ,H) 



2C E 



N " IT { {¥) Vg ~ (r) P B5BMB^/2e + V' A g Al i B H/2e 

+ ( terms independent of N ). (4.12) 



+N e/Cx + a; BC /e - w A (0)/e + 1^1 V c (0) 

§5. Magneto-Coulomb Oscillation in the Ferromagnetic SET 

5.1 The period of the magneto- Coulomb oscillation 

According to the theory of the Coulomb oscillation,!' conductance of an SET in the linear 
response regime is expressed as follows: 

1 1 00 

G= YT Ri +gr £ P cq (N,V g ,H)h(An), (5.1) 
K-qI H T + K T N=Q 

with 

^n = n(Jv,v r g ,fl')-n(Ar-i,y g ,^), (5.2) 

where R l rf mean the tunnel resistances of the left and the right tunnel junctions and P e q(N, V g , H) 
denotes the statistical probability in equilibrium of N electrons being in the island electrode in 
magnetic field H, and is expressed as 

The function h has the form as 

h(x) = r. (5.4) 

1 + exp(x/k B T) 

One can see from the expression of £l(N, V g , H) in eq. (4.12) that, if the value in the square brackets 
in eq. (4.12) is equivalent by a modulus of unity, the sum in eq. (5.1) and, therefore, conductance 
of the SET become equivalent. This fact indicates the periodic nature of the conductance of the 
SET against V g or H. Thus, as a condition which determine the period of the magneto-Coulomb 
oscillation, AH, we obtain the next equation: 

^ h' A 9AmAH/2e - V B g Bt i B AH/2e\ = 1. (5.5) 

As a result, the period of the magneto-Coulomb oscillation is expressed as follows: 

AH= {^) C^ a9a -2CV b9b (5 - 6) 
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5.2 The equi-phase lines in the V g -H diagram 

It is also evident from the expression of Q(N,V g ,H) and eqs. (5.1)-(5.3) that, if the value in 
the square brackets in eq. (4.12) is kept constant with fixed N by tuning both V g and H, the 
conductance of the SET is kept constant. This tuning is realized by the condition such as 

(Ie) Vg ~ (ID V ^ 9BflBH/2e + K9AVBH/2e = constant. (5.7) 
This equation defines the equi-phase lines in the V g -H diagram. Namely they are expressed as 

V g =(^) 2CV ^B-C^V' A g A 

H + constant. (5.8) 

5.3 Expected V g -H diagrams in several ferromagnetic SET's 

Here we discuss some expected V g -H diagrams of ferromagnetic SET's based on the discussions 
in the previous subsections. First we consider the case of Fig. 3(a): the SET is a non-magnetic 
metal/ferromagnetic metal/non-magnetic metal (N/F/N) SET. We assume here the polarization 
V is negative in the ferromagnet. From eq. (5.6) the period of the magneto-Coulomb oscillation is 
AH = (2e 2 / ' hb) I '(Cz\V\g) . The equi-phase lines are expressed as 

J gV'H + constant 

if) fe) 9\ V W H \ + constant 

(if magnetization and H are parallel) 
- (if) (||) g\V\ \H\ + constant 

( if magnetization and H are antiparallel ) . 
For example, suppose one is sweeping a magnetic field from positive strong magnetic field to negative 
one. In positive strong magnetic fields, the direction of the magnetization is parallel to the applied 
magnetic field, V' = +V and the equi-phase lines in this region is expressed by the upper relation 
of eq. (5.9). Once the field is reversed, the magnetization and the field becomes antiparallel as 
long as the magnitude of the field is smaller than the coercive force, H c , of the ferromagnet. In 
this region, V' = —V and the equi-phase lines are expressed by the lower relation of eq. (5.9). At 
H = —H c , the direction of the magnetization changes parallel to the applied field, there occurs an 
instantaneous change in the value of the expression in the braces in eq. (4.12) by 2gV/J,BH c /e, a 
phase jump occurs in the Vg-H diagram, now V' = +V and the equi-phase lines are once again 
expressed by the upper relation in eq. (5.9). The expected equi-phase lines for the peaks (solid 
lines) and valleys (broken lines) of the conductance are schematically depicted in Fig. 3(b). 

Next we examine the case of an F/N/F ferromagnetic SET such as in Fig. 3(c). Here we once 
again assume V is negative in the ferromagnets. In this case, from eq. (5.6), the period of the 
magneto-Coulomb oscillation, AH', becomes AH' = (e 2 / hr) / (C\V\g) , and the equi-phase lines in 
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Fig. 3. Expected V s -H diagrams of two kinds of ferromagnetic SET's. The polarization V of the ferromagnet is 
assumed to be negative. The solid lines are the peaks of the conductance of the SET and broken lines are its valleys. 
H c indicates the coercive force of the ferromagnet. AV g is the period of the Coulomb oscillation. AH and AH' are 
the periods of the magneto-Coulomb oscillations in ferromagnetic SET's of the type (a) and (c), respectively. 



the Vg-H diagram are expressed as 

V g = (jpj gV'H + constant 

-if) (£) 9\V\\H\+ constant 

( if magnetization and H are parallel ) 
V¥){£) 9\V\\H\+ constant 5 ' 10 
( if magnetization and H are antiparallel ) . 

In this case the slope of the lines is in opposite polarity to the case of the N/F/N ferromagnetic 
SET. The expected equi-phase lines in this case for the peaks (solid lines) and valleys (broken lines) 
of the conductance are schematically depicted in Fig. 3(d). 

So far we have assumed that the ferromagnets under consideration have single domains. However, 
in a real ferromagnet of not extremely small size, there exists domain structure in general. In this 
case the magnetization reversal may occur in a finite range of magnetic field, or/and it may occur 
in successive discrete magnetic fields. Partial reversal of magnetization near the junction electrode 
due to the change in the domain structure as well affects the operation of the SET through the same 
mechanism. First, we consider an island electrode that is composed of many magnetic domains. 
We specify the magnetic domain by the number i, name magnetization of the i-th domain mi, its 
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self-capacitance Cj in a magnetic field if.® The shift of chemical potential of the i-th domain, AQ, 
in H is AQi = —VgfiB ^rfffT as cons idered in §2. We assume that, after transfer of electrons among 
the domains to equilibrate their electrochemical potentials in H, excess charge, is accumulated 
in the i-th domain. Then the shift of the electrochemical potential of the i-th domain, Au%, from 
that in zero magnetic field is obtained as Auj( = AQ — eqi/Ci. This quantity should be the same 
among all the domains in equilibrium and represents the shift of the electrochemical potential of 
the island electrode, Auj, in equilibrium. By using a condition, J^f Qi = 0, with N being the total 
number of magnetic domains in the island electrode in H, and noting J^i^i = Cs, we obtain 

Au = _Vm L j : C i rn i .H 

2C S y l m *l 

This means that each domain affects the shift of the electrochemical potential by the factor of Cj/Cs. 
In real devices, we usually observe that capacitances of the junctions are much larger than the other 
contributions to Cy,. In such cases, we can conclude, from the above expression, that the shift of the 
electrochemical potential of the island electrode is almost determined by the magnetizations around 
the junction electrodes because self-capacitances of domains around the junction electrode is much 
larger than the others through the electrostatic coupling with the adjoining counter electrodes. 

Next we examine multi-domain lead electrodes. Here again we specify the domains by the number 
i, name magnetization of the i-th domain m,j and the mutual capacitance between it and the island 
electrode C[. In this discussed in §3, the electrochemical potential of each domain is 

unchangeable by application of magnetic field. The shift of the electrostatic potential of the i-th 
domain, AVi, is expressed as AV{ = —VgfiB ™ \ >m \ ■ The change in the electrostatic potential of 
the island electrode, AV, due to AV^s is expressed as 

2eC s Y \rm\ V ' 

where N' is the number of the domains in the ferromagnetic lead electrodes in H. This expres- 
sion indicates that changes in magnetization around the junction electrodes are most effective in 
changing the electrostatic potential of the island electrode through the large mutual capacitance 
between the island electrode and the domains around the junction electrode. Thus changes in the 
domain structures in the ferromagnetic island and lead electrodes around the junction electrodes 
could cause sequential phase shifts in the V g -H diagram around their coercive forces in sweeping 
the magnetic field. 

On the other hand, if the ferromagnetic component of the ferromagnetic SET is a single crystal, 
the period of the magneto-Coulomb oscillation and the slope of the equi-phase lines in the V^- 
H diagram would show non-monotonic dependences on the magnetic field, which have oscillatory 
structures, in high magnetic fields. It is because the electrochemical potential of the conduction 
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electrons in a metallic single crystal shows an oscillatory dependence on the magnetic field, reflect- 
ing the Landau quantization of the electron orbits in high magnetic fields. 00 Such oscillatory 
electrochemical potential has been actually observed in contact potential measurements. How- 
ever, in a polycrystalline metal such as an evaporated thin film, there exist crystal domains of 
various orientations, and such oscillatory dependence of the electrochemical potential should be 
averaged out. Such effect may be observed in a SET of its island electrode made of a single fine 
particle. 

§6. Consideration on the Chemical-Potential Shift 

As described in the previous sections, magnetic-field-induced shifts of the chemical potentials in 
the ferromagnetic components of the ferromagnetic SET cause the magneto-Coulomb oscillation. 
So far we have considered the chemical-potential shift in a ferromagnet only through the bare 
spin-Zeeman effect. Actually there are contributions to the shift due to other effects. Thus the 
polarization V in the above sections should be replaced by some effective one, V e s- In this section we 
discuss more realistic estimation of the magnitude of the chemical-potential shift in magnetic fields. 
Because numerical calculations based on the one-electron approximation have made a considerable 
success in describing the electronic properties of transition metalsjll) we here discuss it in the 
one-electron approximation. 

6.1 Exchange- correlation correction 

In §2-5, we considered that the bare applied magnetic field alone affects the electron system as a 
perturbation. However, in 3d transition metals, the effective one-electron potential also depends on 
the magnetic field and, as a result, the perturbation is enhanced. Here we examine this effect from 
the point of view of the local spin density approximation (LSDA). As described in §2, population 
of electrons in the states near Fermi energy changes in magnetic fields. As a result, the local 
electron density, n(r), and the local spin density, m(r), change. In the one-electron approximation, 
an electron is assumed to be situated in an effective potential composed of contribution by the 
positive ions and that by the other electrons. Therefore, changes in n and m would make the 
effective potential which an electron feels change. In LSDA© such a change in the effective one- 
electron potential is described by variations in the Hartree potential, vn(r), and the exchange- 
correlation potential, v^ c (r) (a specifying the spin), which are functionals of n and m. Namely, 
an electron feels extra perturbation of Avu(r) + Av° c (r) in magnetic fields. The so-called Stoner 
enhancement of the magnetic susceptibility of the transition metals is considered to be caused by 
this extra perturbation© 

According to the homogeneous enhancement model, which is so far the most successful method 
to calculate the high-field spin susceptibility of the ferromagnetic metals,El'© it is concluded that 
the induced change in n vanishes although that in m is enhanced in homogeneous magnetic fields. 
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This causes an correction to the change in the chemical potential through the variation in v° c {r). 
We obtain the next expression for the change in the chemical potential due to the spin, A( s , in 
this approximation: 

A( S (H) = -V gf i B H(l + T)/2, (6.1) 

where 

T = - - lj (a 3 + a 4 )Xhf • (6.2) 

In the above expression, Xhf and 0^4 respectively represent the high-field spin susceptibility and 
two of the four contributions to the molecular-field coefficient calculated in this approximation in 
ref. 19. Details of the derivation of the above result are presented in Appendix A. 

6. 2 Contribution of the orbital magnetic moment 

Orbital magnetic moments of conduction electrons also contribute to the shift of the chemical 
potential by the change in their magnetic energies in magnetic fields. As it is suggested in the calcu- 
lation of the high-field magnetic susceptibility that the effect of the exchange-correlation correction 
is negligible,© we calculate the orbital contribution to the chemical-potential shift, ziCorb> only by 
the bare perturbation of the external magnetic field. According to the derivation in Appendix B, 
^Corb up to the 1st order in H is given by 

^Corb = ; V j <w\L z \io->, (6.3) 

where i specifies the one-electron eigen state and Co = C{H = 0). As can be seen in the above 
expression, orbital magnetic moments belonging to the states around the Fermi energy cause the 
shift of the chemical potential. Average of them is not necessarily zero in real metals because of 
the spin-orbit coupling. 

Taking into account the results in this and the above subsections, V e R at low temperatures could 
be expressed as 

V cS = V(( )(l+r) + (2/g)< L z > . (6.4) 

Here < L z >o indicates the average of the eigen values of L z of the electron states at the Fermi 
level in zero magnetic field. 

§7. Experimental Observation of the Magneto-Coulomb Oscillation in a Ni/Co/Ni 
SET 

In this section we present experimentally observed results of the magneto-Coulomb oscillation in 
the Ni/Co/Ni SET showing hysteretic behavior against the sweep direction of the magnetic field. 
The device was the same as that used in ref. 4. It is composed of Ni lead electrodes, NiO tunnel 
barriers and a Co island electrode fabricated on a Si wafer with a Ag gate electrode on the rear 
surface of the wafer. The AFM micrograph of the device is presented in Fig. 4. The area of the 
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junctions is about 0.02^im 2 , resistance of the junction, is about 35kf2 and the charging energy, 
e 2 /2Cs, is estimated to be 25-50/ueV. As seen in the figure, there exist many small electrodes of 
the similar shapes with the device itself. However, they are not connected to the measurement 
pads and, therefore, do not contribute to the electric conduction. Further details on the sample are 
described in ref. 4. 



Fig. 4. An AFM micrograph of the Ni/Co/Ni SET device. The electrodes surrounded by a broken-line square 
compose the device. The other small electrodes are dangling and do not contribute to the electric conduction. The 
black spots on both sides of the Co electrode are the tunnel junctions. The gate electrode made of Ag is on the 
rear surface of the chip. 



Measurements were made by sweeping the applied magnetic field slowly and scanning the ap- 
plied gate voltage quickly with the excitation current of 20pA at 12.5Hz to measure the zero-bias 
resistance of the device in a dilution refrigerator at 20mK. The direction of the applied magnetic 
field was parallel to the long axis of the electrodes. 

In Figs. 5(a) and 5(b) are presented the observed results, which show the zero-bias resistance 
of the device by the gray-scale plot on the V g -H diagrams, for both cases of sweeping H to (±) 
directions. (The result in sweeping H to (— ) direction is the same one as in ref. 4.) In both cases 
of sweep directions, the Magneto-Coulomb oscillation is observed clearly. Near zero magnetic field, 
particularly within ±10kOe, there exist many jumps in the equi-phase lines. Moreover, the position 
of these jumps seems to show hysteretic behavior against the sweep direction of the magnetic field. 
Further, there exist solitary phase jumps at high magnetic fields. The magnetic fields at which 
these high-field phase jumps occur are not systematic. They changed sweep by sweep. 

§8. Interpretation of the Experimental Results 

Now we consider the case of the experiment in §7 and in ref. 4 in detail. The key quantity in 
the mechanism discussed above is the spin polarization V of the densities of states. However, V is 
not a directly measurable quantity, and we have to rely on the band calculation. There have been 
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Fig. 5. Hysteretic behavior of the magneto-Coulomb oscillation in a Ni/Co/Ni SET against the sweep direction of 
the magnetic field. Bright (dark) regions correspond to high (low) resistance. The arrow above each frame shows 
the sweep direction of the magnetic field. 

many band calculations on Ni and also several on Co. They have following common features: the 
Fermi energy of the electrons lies just above the majority-spin d bands in both materials, in the 
middle of the minority-spin d bands in Co and near the upper-most of the minority-spin d bands in 
Ni with a much larger density of states (p_) than in Co. Thus it is expected that the polarizations 
are negative in both materials and \Vm\ > | T-'co I - Therefore the expected V g -H diagram is of the 
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type of Fig. 3(d). This is consistent with the experimental observation. 

Some of the band calculations give the densities of states of the majority- and the minority-spin 
electrons at the Fermi energy. Using these values, we calculated Vs in both Co and Ni. The results 
are tabulated in Table I with the values of the densities of states at the Fermi energy per atom 
(p±(°)). In the case of Co the calculated V scatters a little. 



Table I. Densities of states per atom of the majority- and the minority-spin electrons at the Fermi energies and the 
polarizations of Co and Ni. 



material 


p+ (0)* 




V 


Co 


6.64 a 


16.50 a 


-0.43 


(hep) 


2.3 6 


9.0 b 


-0.59 




2.29 c 


10.36 c 


-0.64 


Ni 


2.8 d 


22.3 d 


-0.78 


(fee) 


2.3 b 


22.4 6 


-0.81 




2.860 6 


25.585 e 


-0.80 



* in a unit of states Ryd 1 atom 1 . 

a) from ref. 21, b) from ref. 22, c) from ref. 23, d) from ref. 24, e) from ref. 25. 



On the other hand, sufficient information is not available on the enhancement factor T defined 
in §6 and on the orbital contribution. According to the numerical calculations in refs. 19 and 
20, the factor (03 + a^Xhf m eq. (6.2) is —0.16 for Fe and roughly —0.3 for Ni, which gives 
rather small effect on A£ s . The same conclusion is expected for Co. So we will neglect it in the 
following numerical estimation. The orbital magnetic moment in metals appears due to the spin- 
orbit interaction. However, the magnitude of the spin-orbit interaction in the 3d transition metals 
is not very large; actually, the orbital magnetic moments in these metals are very smallSS'll) So 
we will neglect the orbital contribution as well. 

The period of the magneto-Coulomb oscillation and the slope of the equi-phase lines in the V^- 
H diagram depend also on the parameters in the circuit, namely, the capacitances of the tunnel 
junctions and that of the gate capacitor. In our device, capacitance of the gate capacitor C g is 
estimated to be about 0.60 aF from the period of the Coulomb oscillation. On the other hand, 
the capacitances of the tunnel junctions are not determined accurately in the experiment, because 
the tunnel resistance of the junctions are rather low and smearing of the current (J)-voltage (V) 
characteristic occurs.^ We assumed that the capacitances of both tunnel junctions were equivalent 
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and estimated the magnitude of the sum of the capacitances 2C as 1.6-3.2 fF from the offset voltage 
in the I-V characteristic of the device. Thus we calculate the period of the magneto- Coulomb 
oscillation from the next equation: 



From Table I, Vq — Vm becomes 0.14-0.38. In the experiment, the magnetic field is applied parallel 
to the long axis of the thin film electrodes, so we can ignore the demagnetization effect. Thus the 
expected period of the magneto-Coulomb oscillation AH and the slope &V g /&H become 2.0-11.2T 
and 25-130mV/T, respectively. The observed values, around 2.2T and 120-190mV/T for AH and 
dV g /dH, show a considerable consistency with the theory. 

As observed in Fig. 5 there exist many distinct jumps of the equi-phase lines in both positive 
and negative fields within about ±10kOe. Moreover, they appear mostly on the (+) side when H 
is swept from (— ) to (+) and vice versa. We also observed in ref. 4 a hysteretic tunneling magneto- 
resistance (MR) near H=0Oe: in sweeping H from (— ) to (+), for example, the resistance gradually 
increases from about H=—lkOe till H=+lkOe where the resistance sharply falls. This latter 
phenomenon is ascribed to the so-called magnetic-valve effect: because of the difference of 

the coercive forces of the Co and Ni electrodes, relative orientation of their magnetizations shows 
magnetic-field dependence. This causes magnetic-field-dependent tunneling probability of electrons 
between the two electrodes if the electron spins are conserved during the tunneling processes. From 
the measurement of magnetic induction in ref. 11 of a Ni/NiO/Co tunnel junction, the coercive 
force of Ni is expected to be smaller than that of Co. Therefore the result on the MR in ref. 4 
indicates that, in sweeping H from (— ) to (+), for example, the magnetization of the Ni electrodes 
gradually changes its direction from H = — IkOe and that of the Co electrode changes sharply at 
H=+lkOe in the device. 

Naturally similar hysteretic behavior is expected for the equi-phase lines in the V g -H diagram as 
mentioned in §§5.3. If magnetic field is swept from (— ) to (+), for example, T-^i changes from +Vni 
to —Vm gradually around if=0kOe and that V' Co changes from +Vco to —Vqo rather abruptly 
around H=+YkOe. Moreover, as can be seen in Fig. 4, there is a large Ni electrode on the right 
side adjacent to the narrow Ni electrode. Capacitive coupling between this large Ni electrode and 
the island electrode may cause extra phase jumps or shifts, which are associated with the change in 
the domain structure in this electrode, in a wide magnetic field range. The jumps of the equi-phase 
lines in the V g -H diagram in the region near H=0Oe must be the reflection of these processes of the 
magnetization reversal in both Ni and Co electrodes. On the other hand, the solitary phase jumps 
at high magnetic fields are considered to be ascribed to uncontrollable shifts of the background 
charges. 

As can be seen in the high field range in Fig. 5, the equi-phase lines in the V g -H diagram 




(8.1) 
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show a tendency of leveling off deviating from the straight lines. The origin of this nonlinearity is 
not clear at present. From the semiclassical treatmentlll' 0) of the free energy of a metal in high 
magnetic field (without considering the exchange-correlation correction), there appears a term of 
the 2nd order in H in the non-oscillatory shift of the chemical potential: — (gpBH) 2 (dp/de)/8p, 
where p = p+ + However, based on our self-consistent APW calculation of the density of states, 
this term is estimated to be too small to explain the above mentioned nonlinearity, although the 
energy derivative of the density of states at the Fermi energy in Ni is comparatively very large. 

§9. Possible Applications 

As described in the previous sections the origin of the ferromagnetic SET can be considered 
due to the shift of the chemical potential of the ferromagnet in a magnetic field, and its magnitude 
depends on the electronic properties in the ferromagnet. Taking the discussions in the above sections 
conversely, we see that the ferromagnetic SET can provide a method of measuring a thermodynamic 
coefficient v: 



which is equivalent to another thermodynamic coefficient —(dM/dN)}j,T through the thermody- 
namics. Here M is the magnetization of the ferromagnet. This coefficient contains an information 
on the electronic structure of the ferromagnetic metal complementary to, say, its high-field magnetic 
susceptibility. 

If one uses a SET device tunnel resistance of whose junctions are large, one can determine the 
capacitance of the junction accurately from the current-voltage characteristic of the device as well 
as the capacitance of the gate capacitor, which can be accurately determined from the period of the 
Coulomb oscillation. Therefore accurate measurement of v is possible. The measurement procedure 
is as follows: the appropriate device is a symmetric N/F/N ferromagnetic SET of the type depicted 
in Fig. 3(a) of rather high tunnel-junction resistance. After applying the SET a strong magnetic 
field, for example, in (+) direction in order to align the magnetization of the island electrode to 
the (+) direction, one decreases the applied magnetic field to near zero field and measure the 
slope (dVg/dH) of the equi-phase line in the Vg-H diagram. By using the numerical values of the 
capacitances C and C g , which are determined in other measurements, v is determined from the 
slope of the equi-phase line by the next relation: 




(9.1) 





Within the approximation in §6, where the exchange-correlation enhancement of the perturbation 
is included only in the spin part, the expression of v becomes as 



v = -VgpB (1 +T)/2-p B < L z > . 



(9.3) 
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Apparently T contains an information on the electronic structure complementary to the Stoner 
enhancement factor. This method would provide a new tool to investigate electronic structure 
in ferromagnets complementary to the methods such as the field emission, the spin-resolved pho- 
toemission, the magneto-optic Kerr-effect, the spin-polarized tunneling^!* and the high-field sus- 
ceptibility measurement. In the present stage of experimental technology, fabrication of N/F/N 
ferromagnetic SET's of variety of ferromagnets with appropriate device parameters is not difficult. 
Actually, as a first trial, we have made an Al/Co/Al SET with the tunnel barrier of AI2O3 and 
estimated the magnitude of v of Co from the measurement in the magnetic-field range where Al is 
normal-conducting@ 

§10. Conclusion 

We have shown that magnetic-field-induced variations in the chemical potentials in ferromagnetic 
components of the ferromagnetic SET changes the free energy of its island electrode and cause a 
conductance oscillation (magneto-Coulomb oscillation) similar to the Coulomb oscillation by sweep- 
ing an applied magnetic field. The recently observed magneto-Coulomb oscillation in a Ni/Co/Ni 
ferromagnetic SET can be consistently interpreted in this mechanism. Some related phenomena, 
magnitude of the magnetic-field-induced variation of the chemical potential of the conduction elec- 
trons in a ferromagnetic metal, and an application of the ferromagnetic SET as a tool to investigate 
the electronic structure in ferromagnetic metals have also been discussed. 
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Appendix A: Calculation of the Exchange-Correlation-Correction Factor 

In LSDA the total energy of the system without external magnetic field is written asE^ 

E (n,m) = -Ze 2 J2 I ^ r T~Z^R\ + T / ^ r I dV ^j^ + T s (n,m) + E xc (n,m), (A-l) 

it 

where n(r) is the electron density. The first term denotes the Coulomb energy due to the nucleus 
with charge Ze at the site R and the second term denotes the Hartree energy. The third and 
fourth terms denote the kinetic energy and the exchange-correlation energy, respectively, which are 
functionals of the electron and spin densities, n and m. The exchange-correlation energy £ 1 xc (?i, 771) 
is assumed here to be known. From eq. (A-l) the unperturbed Hamiltonian for an electron with 
spin a is obtained:© 

h 2 

H° a = -—V 2 + v(r) + VH(r) + v° c (r), (A-2) 



20 



where 

v{r) = -Ze^ ¥ L ^y (A-3) 

^(r) = e 2 / dV^l, (A-4) 
J |r — r | 

r -N _ SE xc (n,m) 5E xc (n,m) 

Here ft denotes the Planck constant and m e means the electronic mass. In order to make it 
more complete, a term representing the spin-orbit coupling operator should be added to the above 
expression of H®. In an external magnetic field H, represented by the vector potential A, the 
Hamiltonian may be written in the form: 

n (7 = n° a + -^A 2 + n' CT . (A-6) 

2m e c z 

Here TC' a represents the enhanced perturbation and is given by 

K = -Wff(r)/2 + Av n {r) + Av° c (r) - —A • V. (A- 7) 

m e c 

In this expression, the first term represents the bare spin-Zeeman effect. The second and third 
terms arise since magnetic field changes the local electron and spin densities n(r) and m{r). The 
last term represents the orbital term, the effect of which is discussed in Appendix B and neglected 
here. The induced changes in the Hartree and the exchange-correlation potentials make energies 
of electrons shift from those including the bare spin-Zeeman effect, and cause a correction for the 
chemical potential shift. 
The Helmholtz free energy of the conduction electron system is given by 

F = N( + Tr$(H), (A-8) 

where iV is the number of conduction electrons in the metal, £ is the chemical potential and 

$(H) = -k B Tln[l+exp(( -H)/k B T}. (A-9) 

The chemical potential is determined by the next condition so as to minimize the free energy: 

dF 

af-o. (a.io) 

According to the Peierls' technique, ^(TC) may be expanded as follows ,13 §1) 

Tr $(H) = Tr <2>(H ) + £ <2>'(e A ) < X\K\X >+•••, (A-ll) 

A 

where the representation in which Tt® + e 2 A 2 /(2m e c 2 ) is diagonal is used. Based on this expan- 
sion we calculated the shift of the chemical potential between with and without H! a . Discarding 
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negligibly small terms, the result up to the 1st order in H is as follows: 

Co) • I ■ . 




(A-12) 

< icr|Z\u H + Av° c \ia > . (A-13) 

The second equality holds because, up to the 1st order in H, the eigen states of 7i° a can be used in 
the sum. At low temperatures the expression becomes as 

A( s = -V gf i B H/2 + - d/(< T ~ < »^VH + 4< c |i<i > . (A-14) 

The first term denotes the contribution by the bare spin-Zeeman effect. The second term indicates 
that average energy shift due to Avu(r) and Av° c (r) for the electron states around the original 
Fermi energy composes the correction term. The same expression could be obtained at high mag- 
netic field and at low temperature based on the semiclassical theory of the de-Haas van Alphen 
effect© 

The changes in n{r) and m(r), which cause Av-ft(r) and Av° c [v), are induced by the shifts of 
the electron energies, and the electron energies are determined by the Hamiltonian which includes 
Avn(r) and Av° c (r). Therefore, Avn(r) and Av° c (r) should be calculated in a self-consistent 
way. In the homogeneous enhancement model,© induced changes in m(r) and n(r) by an external 
magnetic field with a wave-number vector, q, are assumed to be such as Am(r) = aqAm s (r) 
and An{r) = 0qAn s (r), where aq and (3q are q-dependent constants, and Am s and An s are the 
unenhanced spin-density and charge-density shifts, respectively. Namely, the induced changes in n 
and m are assumed homogeneously enhanced independent of the position r. The authors of ref. 
19 solved the 1st order perturbation equations for Am(r) and An(r) self-consistently at finite q in 
LSDA. Then they calculated ao, the enhancement factor for the change in m(r) in a homogeneous 
magnetic field, by taking the limit of q— > to obtain the high-field spin susceptibility XhS- By the 
similar calculation, /3q, the enhancement factor for the change in n{r) in a homogeneous magnetic 
field, can be easily obtained: 

/3o = 0. (A- 15) 

This means that, in a homogeneous magnetic field, local electron density, n(r), does not change. 
Because the Hartree potential is a functional of n{r) alone, Avn(r) vanishes and one should consider 
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only Av° c (r) in TL' a . Moreover, because An(r) = 0, Av% c (r) is simplified in this case as 

Av° c {r) = J d 3 r(L(r) + aJ(r))Am{r), (A-16) 

where L(r) and J(r) are the second functional derivatives of the exchange-correlation energy 
E xc (n,m) as defined in ref. 19. Substituting (A-16) in (A-14) and using a an d Xhf obtained 
in ref. 19, expressions for the enhanced (suppressed) shift of the chemical potential (6.1) and (6.2) 
are obtained. 

Appendix B: Calculation of the Orbital Contribution 

Here /\Corhn the contribution of electron orbit to the shift of the chemical potential in magnetic 
fields, due to the bare perturbation of the external magnetic field is calculated. The Hamiltonian 
representing the perturbation related to the orbit may be written in the form: 

*U~^.V CB-l) 
= VbHL z . (B-2) 

Here the magnetic field is assumed in z direction and the usual gauge A= (Hxr)/2 is used. This 
term adds the next term to Tr$(7Y) in (A-ll): 

£ <P'(e A ) < A|HU|A > ■ (B-3) 

A 

As a consequence, it causes a shift of the chemical potential of the next expression, which is additive 
with A( s : 

A( mh = -^E d/(, T C0) < A "^ A > • (B-4) 

P++P-X d£ A 

Considering up to the 1st order in H, the eigen states of can be used in the above expression. 
Thus, we obtain eq. (6.3). 
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